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ABSTRACT 


a ki,, ,vf th* Qn arp Station reauires numerous construction flights by the Space 
Shuttte A particularly challenging problem is that of control of each intermediate station 
confieurati^ OrtSter is approaching it to deliver the next component 

The necessary braking maneuvers cause Orbiter thruster plumes to impinge on 

arrays. This in tum causes both overall attuude errors and 

excitation of flexible-body vibration modes. These plume loads are predicted to lead 
CMo irat on XnTthe approach of the Orbiter to the SC-5 stauon configuration 
SSt use of the station RCS jets for desaluration. They are also expected to 
lead to significant excitation of solar array vibrations. 


It is therefore of great practical importance to investigate the effects of plume loads on 

dynamics i e' those which are strongly excited by the control inputs or plume disturbance 
forces and which strongly couple with the measured outputs. The particular technique to 
£ 'SJta ^hfs study is the subsystem balancing approach which was previously 
Severed by ^ the present investigator. This method is very efficient computationally 
FurUiennore it gives accurate results even for the difficult case where the structure has 
many closed-spaced natural frequencies, when standard modal truncation can give 
misleading resuUs. Station configuration SC-5 is a good example of such a structure. 
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INTRODUCTION 


Assembly of the Space Station requires numerous construction flights by the Space 
Shuttle, regardless of the details of the station configuration that is finally selected. A 
particularly challenging problem is that of control of each intermediate station 
configuration when the Shuttle Orbiter is approaching it to deliver the next component. 
The necessary braking maneuvers cause Orbiter thruster plumes to impinge on the 
station, especially its solar arrays. This in turn causes both overall attitude errors and 
excitation of flexible-body vibration modes. These plume loads are predicted to lead to 
CMG saturation during the approach of the Orbiter to the SC-5 Space Station Freedom 
(SSF) configuration, necessitating the use of the station RCS jets for desaturation. They 
are also expected to lead to significant excitation of solar array vibrations. 

It is therefore of great practical importance to investigate the effects of plume loads on the 
flexible dynamics of station configuration SC-5 as accurately as possible. However, this 
system possesses a great many flexible modes (89 below 5 rad/s), making analysis time- 
consuming and complicated. Model reduction techniques can be used to overcome this 
problem, reducing the system model to one which retains only the significant dynamics, 
i.e. those which are strongly excited by the control inputs or plume disturbance forces and 
which strongly couple with the measured outputs. A considerable amount of work has 
been carried out on the topic of model reduction in the past. Well-known methods 
include modal truncation [1], based either on the natural frequencies of the stmcture or its 
modal costs, and balancing [2] of the entire structure and then truncation to retain a 
dominant model for it. An advantage of the balancing approach is that it typically yields 
a more accurate reduced-order model than does simple modal truncation. This is 
particularly true when the structure possesses clustered natural frequencies, as is often the 
case for realistic flexible space structures. However, the disadvantages of balancing are 
its high computational cost, possible numerical sensitivity problems resulting from the 
large matrices being operated on, and the difficulty involved in providing a physical 
interpretation for the resulting balanced "modes". 

The purpose of the work reported here is to investigate the performance of the alternative 
subsystem balancing technique when used to study the plume impingement problem on 
SC-5. This method, introduced in [3][4], further developed in [5], and then applied to a 
simplified SSF model in [12], retains the desirable properties of standard balancing while 
overcoming the three difficulties listed above. This is achieved by first decomposing the 
structural model into subsystems of highly correlated modes, based on the modal 
correlation coefficients derived in [4] from the controllability and observability 
Grammian matrices [6] of the structure. Each subsystem is approximately independent of 
all others, so balancing each separately and concatenating the dominant reduced-order 
models obtained yields roughly the same result as balancing the entire structure directly. 
The computational cost reduction produced by this subsystem technique is considerable: 

an operation count reduction by a factor of roughly r if the system decomposes into r 
equal subsystems. The numerical accuracy of the resulting reduced-order model is also 
improved considerably, as the matrices being operated on are of reduced dimension; this 
avoids the numerical conditioning problems noted in [8] [9] for standard balancing. 
Furthermore, the modes of the reduced model do now permit a clear physical 
interpretation. This is a consequence of the fact that each correlated subsystem must 
necessarily only include modes with close natural frequencies. The balanced modes of 
each subsystem are therefore, to first order, linear combinations of repeated-frequency 
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modes, and so can themselves be taken as an equally valid set of physical modes. 
Balancing the entire structure, on the other hand, combines modes of widely differing 
frequencies, making interpretation difficult. 

The graphs and tabulated results given in this report demonstrate the improvements 
achievable by using the subsystem balancing technique for model reduction. In fact, the 
original 94-mode model for the SC-5 SSF configuration was reducible to a 10-mode 
model at the cost of a modeling error of less than 2%. This then allows analysis and 
simulation of the effects of plume loading to be carried out much more efficiently than if 
the entire model were used. 

THEORETICAL BACKGROUND 

Consider an n-mode model for the structural dynamics of a modally damped, non- 
gyroscopic, non-circulatory FSS with m actuators and p sensors, not necessarily 
collocated. This model can be written in modal form [1] as 

r\ + diag(2^co, )T\ + diag(Q) f )t| = fiu, ^ 

y = c r i\ + c d x\, 

where t\ is the vector of modal coordinates, u that of applied actuator inputs and y that of 
sensor outputs, and (0 { and £, are the natural frequency and damping ratio of the i th 
mode, respectively. For the typical FSS [7], the {£", } are quite low (e.g. 0.5 %), and the 
occur in clusters of repeated, or nearly repeated, frequencies as a result of structural 
symmetry. 

Defining the state vector x = (f} t , ( 0 ^ ,•• co„rf n ) T for this structure yields the state 

space representation x = Ax + Bu, y = Cx, where A = blkdiag(A)j, B — (B ] ,*••, B n ) and 
C = (Cp-'-.C,), with 


A ; = 


(- 2 ^ -G>A = fb 


0 ): 


0 


q j and C, =(c n c */©,); 


( 2 ) 


bj is the i 1 * 1 row of B , and c n and c^, are the i^ 1 columns of C r and C d , respectively. 
The problem studied here is that of obtaining a reduced-order model 


= 4 X r + B r U. (3) 

y r = c r \ r 

for this structure for which the normalized output error 


$ 2 = JIy(0-y r (0l2^/J||y(0||2<* 
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is aeecptohiy small. Of course, the size of 8 will depend on the order, n r , chosen for the 
reduced model. A good model reduction procedure should ideally provide information 
allowing an intelligent choice for n r to be made so as to achieve a specified 8 value. 

Two techniques for model reduction that have been extensively studied are those of 
modal truncauon and internal balancing. The new method implemented in this report 
subsystem balancing can be regarded as a generalization of the two established 
techniques; it includes both as special cases, but permits a wide range of other solutions 
also It is this additional freedom that allows subsystem balancing to produce superior 
results to those obtainable by means of the traditional methods. F 

Model reduction by subsystem balancing proceeds by first dividing the given structure 
into subsystems of highly correlated modes. The dimensions of these subsystems are 
determined by the choice of a correlation coefficient threshold value by the analyst. This 
quantity is used as follows. Two modes which have nearly equal natural frequencies and 
non-orthogonal mode shapes will be approximately in resonance, and so can strongly 
excite each other An input which excites the first mode will thus typically also indirectly 
excite the second. It is therefore important to treat the two modes as a single unit, or 
subsystem, when considering whether or not to retain them in a reduced-order model. 
i>uch a pair of close modes can be shown to have a high correlation coefficient; the basis 
ot subsystem balancing is therefore to define subsystems made up of all modes which 
have correlation coefficients greater than the defined threshold value. It can be seen that 
selecting a threshold of 0 implies that all modes are deemed to be correlated, i.e there 
will be a single subsystem" which is in fact the entire system. This is equivalent to the 
case of standard balancing. Conversely, choosing a threshold of 1 leads to all modes 
being deemed to be uncorrelated, and so dealt with in isolation; this amounts to the modal 
truncauon technique. 


Each of the subsystems obtained for the chosen threshold value is then balanced 
independently, and a reduced-order model for it generated by deleting all balanced states 
corresponding to Hankel singular values [2] below some specified threshold. Note that 
the singular value weighting described in [10] could be applied, if desired, without 
changing the argument in any way. Similarly, frequency weighting of the Hankel 
singular values can easily be incorporated to deal with input signals which have a known 
frequency spectrum. This is actually done in the present application, where the inputs are 
steps (representing thruster firings) rather than the impulses classically considered in 
model reduction problems. The resulting reduced-order subsystem models so obtained 
me then combined to yield a dominant, approximately balanced, reduced-order model for 
the full system. 


USER INTERFACE 


This section describes the user interface to the model reduction package which was 
developed as part of this contract. This software consists of a library of Matlab m- 
functions, with mrmain3 calling all the other functions internally. The package is installed 
on the Sun SparcStation 2 carmel in the Interactive Analysis Laboratory in Building 16 
and has also been provided in a Macintosh version to personnel at JSC and Dynacs 
Engineering Co. The information that follows details the user interface for mrmain3. This 
function, together with the purpose-built second level functions that it calls, have 
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extensive in-line documentation, facilitating future use and/or modification. Listings of 
all these functions are provided in the report [12], 

Input arguments 

om: The natural frequencies (rad/s) of the structure, input as either a row or column 
vector. Any rigid-body modes must precede the flexible modes and be represented by 
hard zero frequencies. 

phia: The influence matrix, in mass-normalized coordinates, corresponding to the 
specified actuator locations. If the structure has n modes and m actuators, phia will be an 
(n x m) matrix. 

phis: Similar to phia, but for sensor stations or positions of outputs of interest (e.g. solar 
array tips). 

Prompts to user 

The following prompts are generated by mrmain3 when running. They allow each run to 
be customized as desired by the user. 

Output the time taken for each step?: The time required for each matrix decomposition, 
etc., is output to the screen if requested. This allows the progress of the model reduction 
procedure to be monitored, as well as giving an indication of which steps are the most 
computationally intensive. 

Vectorize? ( Faster , but requires more storage): In Matlab, for loops are typically an 
order of magnitude slower to execute than the equivalent "vectorized" operation. For 
instance, s=0; for i=1:n, s=s+x(i); end; runs considerably slower than does s=x*ones(n,1). If 
vectorization is requested, computation of the system Grammian matrices and correlation 
coefficients is put into the form of vector-matrix operations rather than loops; this is 
indeed considerably faster, but requires some additional temporary storage arrays. 

Structural damping ratio, % ( default is 0.5%): The specified damping ratio is applied 
uniformly to all flexible modes of the full structural model. 

Print frequencies in Hz?: The mean frequency of each subsystem can be output in either 
rad/s or Hz, as desired. 

Desired controllability threshold?: This threshold value is used to determine which 
modes are correlated in a controllability sense. The system is then broken down into 
disjoint sets of modes (subsystems), where modes with a controllability correlation 
coefficient greater than the specified threshold are deemed to be correlated. Taking a 
threshold value of 0 implies that all modes are considered correlated, i.e. the method 
reduces to standard balancing. Conversely, a threshold value of 1 implies that no modes 
are taken together: this is modal truncation. Intermediate values allow the dimensions of 
the resulting subsystems to be specified to a large extent; reducing the threshold reduces 
the number of subsystems, so increasing their dimension. The program allows the user to 
try various different values for this threshold until he is satisfied with the resulting set of 
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subsystems. At this point, a negative value is input to cause the program to proceed to the 
next step. 

Desired overall threshold?: This threshold is used in a similar fashion to the 

controllability threshold, but both controllability! and observability are now taken into 
account. This yields the final subsystem distribution output by the program (in modemap) 
and used to obtain the reduced-order model. 

Compare step responses?: If requested, the step responses of the full and reduced-order 
models are computed, plotted, and the relative differences (i.e. reduced-order model error) 
output for each input-output channel. 

Compare frequency responses?: If requested, the Bode amplitude plots of the full and 
reduced-order models are similarly computed and plotted, and the relative differences (i.e. 
reduced-order model error) output for each input-output channel. 

Desired truncation measure?: Two types of measure can be used to define the number of 
modes retained in the reduced-order model. If a positive integer is input, this is taken to 
be precisely the desired reduced model order. On the other hand, if a real number in the 
interval [0, 1 ) is input, this is taken to be the desired relative error in the reduced-order 
model step response, and the model order required to achieve this is computed. (Note 
that this second option is only an approximation, and should be treated as such.) The 
program loops over this step, allowing the user to try a sequence of various different 
reduced model orders to determine which one is most satisfactory for his purposes. This 
is quite a rapid operation when using mrmain3, as nearly all computations need only be 
done once. Essentially the only step that must be repeated when trying a new model 
order is that of calculating the step and/or Bode plots of the reduced-order model. 

Output arguments 

am, bm, cm: The reduced-order state-space model matrices {A r , B r , C r } obtained. 

modemap: This matrix specifies which physical modes are grouped into which 

subsystems in the decomposition based on overall correlation coefficients. In particular, 
the i 1 * 1 column of modemap lists the modes making up the i^ 1 of these subsystems. 

DISCUSSION OF RESULTS 

Results will now be provided which illustrate the behavior of the subsystem balancing 
technique when applied to a structural model [11] of the SC-5 configuration of Space 
Station Freedom. This structure possesses light damping (estimated to be 0.5% of 
critical), and a large number of closely-spaced vibration modes (94 flexible modes in the 
model considered, of which 89 are below 5 rad/s). The model has 12 inputs: 6 Reaction 
Control System (RCS) thrusters, and 6 disturbance inputs which represent the lumped 
forces and torques about all three axes that result from Shuttle Orbiter plume 
impingement during approach. Note that the Control Moment Gyroscopes (CMGs) that 
would also be used on SC-5 do not need to be considered as inputs here: this is because 
CMGs incorporate torque shaping techniques so as to not significantly excite flexible 
vibration modes. The measured outputs are the 3 angular rates sensed by the rate gyros 
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on the station. (The movements at other positions of interest, for instance the solar array 
tips, could also be considered if desired; the method remains exactly the same.) 

The three numerical parameters which the user must select when running mrmain3 are: the 
controllability coefficient threshold value; the overall (controllability and observability 
combined) coefficient threshold; and the order of the final reduced-order model. It 
should be noted that, as previously pointed out, the program allows the user to try a 
sequence of values for each of these quantities until satisfactory results are obtained. The 
effects of different choices for these three parameters will now be examined. 

Figure 1 illustrates the role of the controllability threshold coefficient p cl in determining 
the subsystem decomposition. The two solid graphs give the maximum and minimum 
subsystem dimensions that result for values of p ct ranging between 0 (standard balancing) 
and 1 (modal truncation). The dashed graph then shows the number of subsystems 
obtained for each threshold value. It can be seen that the system decomposition does 
indeed change, as expected, from that of standard balancing (the only subsystem is the 
entire model, with 94 modes) to that of modal truncation (there are 94 single-mode 
subsystems) as the threshold increases from 0 to 1. It can be noted that the evolution of 
subsystem dimensions is somewhat discontinuous: for instance, large changes occur for 
thresholds between 0.05 and 0.15, whereas there are hardly any differences between 0 55 
and 0.70. A consequence of this is that it is not always possible to find a threshold value 
which will yield a particular maximum subsystem order. However, it is always possible 
to obtain a good working value which gives a totally acceptable subsystem partition. 
From practical experience, it has been observed that a good choice of threshold value is 
generally one which gives the maximal subsystem dimension approximately equal to the 
number of subsystems. The reason for this is that it provides a good balance between 
having subsystems which are too large (and so susceptible to numerical problems) and 
ones which are too small (and so neglect appreciable cross-mode interaction). For the 
system studied here, such a choice can be seen to be p ct = 0.2, giving 27 individual 
subsystems (some of which consist of single modes), and a maximum subsystem 
dimension of 17. 

Figures 2-4 similarly illustrate the role of the overall correlation coefficient threshold p ol . 
These figures can be seen to be broadly similar to Figure 1. The main difference is that 
the choice made for p cl affects the initial, controllability-based, subsystem 
decomposition, which is then used to compute the final subsystem decomposition based 
on p ot . Figures 2-4 demonstrate this by letting p ot run through all possible values 

between 0 and 1 for p cl equal to 0.1, 0.2 or 0.3, respectively. The differences between the 
three plots can clearly be seen: the main distinction is that the low controllability 
threshold value of 0. 1 tends to give rise to larger final subsystems. This makes sense, as a 
lower threshold leads to more modes being considered sufficiently correlated that they 
must be grouped in the same subsystem. A rule-of-thumb for overall threshold that has 

generally been found to work well is simply to take Po, = p ct . 

Table 1 then quantifies the effect of different choices for the reduced model order n r . In 
all cases, the controllability and overall thresholds are taken to be equal, p ol = p cl — p t , 
and this value is allowed to range over values between 0 (balancing) and 1 (modal 
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truncation). The numbers that are tabulated are the step response relative error norms, as 
defined by (4), for the resulting reduced-order models. It is significant to note that the 
new subsystem balancing technique gives better results than those of standard balancing 
or modal truncation for all the reduced orders considered. For instance, reducing the 
original 84-mode system to a 10-mode model by subsystem balancing with p, = 0.2 gives 
rise to a relative error of less than 2%. By contrast, the error obtained for a model of the 
same dimension by means of balancing is over twice as large, and that produced by modal 
truncation is nearly three times as great. 


TABLE 1.- STEP RESPONSE RELATIVE ERRORS 
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Figure 5 illustrates this further by plotting the step responses of both the full system and 
the 10-mode reduced-order model obtained for p, — 0.2. It can be seen that the plots are 
nearly indistinguishable, despite the nearly order-of-magnitude reduction of system 
dimension. Figure 6 confirms this by plotting the step response error between the 
reduced model and the original system: note the different y-axis scale factors. For 
comparison. Figure 7 gives the error that would be obtained for a 10-mode model 
generated using modal truncation; again note the change in y-axis scale. The increase in 
this error over that obtained by subsystem balancing is clear. Finally, Figure 8 overlays 
the frequency response plots of the full system and the 10-mode reduced model obtained 

by subsystem balancing for p, = 0.2. (These plots are weighted by a 0)~' term, so as to 
reflect die fact that a step input is being applied.) The close match between the two 
models is clearly visible in the frequency domain, just as it was in the corresponding time 
histories. It is also easy to identify which flexible modes are strongly excited by a 
prolonged, step-like, plume excitation force. 
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Max and min size and number of subsystems Max and min size and number of subsystems 


Subsystem balancing on SC- 5 



Controllability correlation threshold 


Figure 1.- Effect of varying controllability threshold level 


Subsystem balancing on SC-5, ro-c = 0.1 



Overall correlation threshold 


Figure 2.- Effect of varying overall threshold level for p cl = 0.1 


EKE 


3 ‘ 


Subsystem balancing on SC-5, ro-c = 0,2 



Figure 3.- Effect of varying overall threshold level for p cl = 0.2. 



Figure 4.- Effect of varying overall threshold level for p ct = 0.3. 
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xlO- 6 Step responses, nrom = 1 0; input 1 , output 1 



Figure 5.- Step responses of full and reduced models. 

xlO** Step error, nrom = 10; input 1, output 1 



Figure 6.- Step response error of reduced model. 


34-12 




xlO 7 Step error, nrom = 10; input 1, output 1 



Figure 7.- Step response error of modal truncation model. 


Frequency responses, nrom = 10; input 1, output 1 



Figure 8.- Frequency responses of full and reduced models. 
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